Atomic decoration for improving the efficiency of field electron emission of carbon nanotubes 
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The field electron emission from the single- walled carbon nanotubes with their open ends terminated by -BH, 
-NH, and -0 has been simulated. We find that -BH and -NH suppress the apex-vacuum barrier significantly 
and lead to higher emission current in contrast to the -O terminated structure. The calculated binding energy 
implies that the carbon nanotubes terminated with -BH and -NH are more stable than those saturated by oxygen 
atoms or by hydrogen atoms. The simulation shows that the most probable orientation of the emission beam has 
correlation to the atomic structure of the apex and is field-dependent. 
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I. INTRODUCTION 

The field electron emission (FE) from carbon nanotubes 
(CNTs) has found its applications in flat panel displays dl] 
and in molecular sensors. (2, S 0] It has also shown poten- 
tial applications in miniature highbrightness electron sources 
for electron microscope [5] and in parallel e-beam lithogra- 
phy system. (6D One of the central problems in these appli- 
cations is to optimize the structure to achieve stronger FE in 
lower macroscopic field. It was once believed that the large 
aspect ratio of a CNT that would lead to large local field en- 
hancement (LFE) at the apex of the CNT was the major reason 
responsible for the superior FE, hence it would be straight for- 
ward to improve the FE efficiency by increasing the length of 
the CNT. But this common view has not been verified satis- 
factorily by either experimental or theoretical studies. Recent 
simulations have shown that the charge accumulation in both 
the apex and the body of the CNT can significantly modify 
the apex- vacuum barrier (AVB) and thereby the FE charac- 
teristics. |0 [sh A recent simulation found that the field en- 
hancement factor for the open single- walled CNT (SWCNT) 
is much smaller than the expected value of the metal rod 
model; on the other hand, the change of the AVB has more 
pronounced effect to FE. Q2U 

Since a CNT emitter has only decades of atoms at the apex, 
the atomic structure of the apex would have strong influence 
on the FE process. The following effects have been observed 
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in experiments: the hydrogenation of the tube wall trans- 
forms a metallic CNT to a semiconductive one, II 1011 O2 ex- 
posure increases the turn-on field of SWCNTs and decreases 
the FE efficiency, [11] and the adsorption of H2O enhances 
the field emission current. [12] However, experimental ob- 
servations for the adsorbate effects so far have not been con- 
clusive, fill IT3I |3 To understand the dependence of FE 
upon the atomic structure of the apex, careful simulations via 
the density functional method (DFT) have been carried out. 
11151 lid 1 1 VL ll 8l Il9l 1 2 (A |2 ill However, there are contradictory 
conclusions about the effect of adsorbates. Zhou et al. II22I1 
and Kim et al. |0 HH] obtained the local density of states 
(LDOS) at the apex by the ab initio methods and found that 
the LDOS at the charge-neutrality level is suppressed by the 
hydrogen. They therefore concluded that hydrogen adsorption 
reduces the FE current density. By contrast, Mayer et al. cal- 
culated the AVB using a dipole and point charge model. [0] 
They acclaimed that the apex- vacuum barrier is reduced by the 
presence of the hydrogen, and thereby concluded that hydro- 
gen adsorption would enhance the FE current density. Mayer 
recently improved the model and illustrated the electrostatic 
potential around the carbon nanotube. 1 10] More careful stud- 
ies on this topic would obviously be useful. 

Only recently, it is possible to tackle the SWCNT of realis- 
tic size in the FE conditions by a multi- scale method involv- 
ing quantum mechanics and molecular mechanics, |24|] In 
the present paper, we have adopted this method to simulate 
the FE from the SWCNT of realistic length (one microme- 
ter) with different atomic decorations at the open end of the 
tube. Obviously the AVB and thereby the FE characteristic 
will be strongly affected by the electron transfer between the 
carbon atoms and the adsorbent atoms. The atomic decoration 
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TABLE I: Electronegativity of related atoms in Allen electronegativ- 
ity scale. 



atom 


H 


B 


C 


N 





X s 


2.300 


2.051 


2.544 


3.066 


3.610 



in the apex will induce dipoles due to the geometric symme- 
try breaking in the axial direction. If the dipole has its pos- 
itive end outward to the vacuum (positive dipole), it tends to 
suppress the AVB; otherwise (negative dipole), it inclines to 
raise the barrier. This simple argument suggests that the car- 
bon dangling bonds in the open end of the SWCNT should be 
saturated by atom with electronegativity (Xs) lower than the 
carbon atom. For instance, the hydrogen terminated SWCNT 
has a positive dipole because hydrogen has lower Xs than car- 
bon. [9] The electronegativities of hydrogen, boron, carbon, 
nitrogen, and oxygen in Allen electronegativity scale are given 
in TABLE H 

At room temperature, the electrons in the vicinity of the 
apex are most relevant to the emission. When the apex has 
higher density of electrons, there are more incident electrons 
hitting on the AVB and a stronger emission current could be 
expected. To attract electrons from the tube column to the 
apex, it would be useful to saturate the CNT first by atoms of 
higher X s then followed by atoms of lower X s . Another con- 
sideration is the structural stability that is one of the most con- 
cerned properties in the applications. In this paper, we should 
consider the diatom ions -BH and -NH as adsorbates and com- 
pare them with the oxygen ion. Our simulation should show 
that the adsorbates of -BH and -NH have larger binding energy 
than both oxygen and hydrogen. 

We should adopt the multi-scale method (MHHl to simulate 
the (5, 5) SWCNT with its open end saturated in three kind 
of ions: -BH, or -NH, or -O. In Section II, the method of the 
simulation is reviewed briefly. The charge density and elec- 
trostatic potential in the absence of applied field is given in 
Section III. In Section IV, the field-depending AVB is illus- 
trated. The most probable emission path is studied and the 
emission characteristics of the (5, 5) SWCNTs saturated with 
-BH, -NH and -O are presented in Section V The last section 
is devoted to the discussions and conclusions. 



II. SIMULATION METHOD AND THE ENDING 
STRUCTURES 

The concept of multi- scale coupling is important for the 
simulation of huge systems that are sensitive to all scales of 
the systems. ||25l |26|] The CNTs for the purpose of FE are 
typical multi-scale systems. A simplified schematic setup of 
the FE system is shown in FIG. [T] in which two black plates 
are the cathode (left) and anode (right); a CNT is mounted 
vertically on the cathode. When a voltage is applied to the 
two metal plates, the electrons have opportunity to emit into 
vacuum through the apex of the CNT by quantum tunneling. 

In experiments, the length of CNTs is usually in microme- 
ters, while the radius is in nanometers. For instance, the (5, 5) 
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FIG. 1: The schematic setup of field electron emission. The arrow 
lines represent the applied electric field. 

type SWCNT of one micrometer long consists of about 10000 
carbon atoms. Limited by the computational efficiencies and 
resources, all ab initio studies so far can only simulate the 
local properties involving hundreds of carbon atoms. As the 
electronic properties are sensitive to both the detailed atomic 
arrangement (i.e., the location of defects and adsorbates) and 
the electron states, which would extend through the whole 
tube, it is a big challenge to simulate the whole tube that has 
a length in micrometer scale. 

Since electrons are emitted from the apex of the SWCNT by 
quantum tunneling, the apex part must be treated by quantum 
mechanics. The part on the substrate side mainly affects the 
FE through Coulomb potential of the excess charge, so that 
it can be treated by a semi-classical method. $H [24]] There- 
fore we divide the SWCNT into a quantum region and a semi- 
classical region. The quantum region is simulated on atomic 
scale where the density matrix of electron is obtained quan- 
tum mechanically. In the semi-classical region, the Coulomb 
potential is governed by the Poisson's equation. 

It should be noted that even in the semi-classical region the 
electron energy band structure originated from quantum me- 
chanics should be taken into account. The excess charge dis- 
tribution in the semi-classical region has simple solution only 
for the simple band structure. For the (5, 5) SWCNT, there 
are both experimental and theoretical evidences for the con- 
stant density of states (DOS) in the vicinity of neutrality level. 
J27L Hsl [23, [5(1 Hence the semi-classical region should be 
limited in the region where the constant DOS is valid. To our 
experience, this means that the length of the quantum region 
should be over 100 nm. 

In the present simulation, the quantum region extends 123 
nm and contains 10000 atoms. It is still too big for an ab initio 
simulation. The quantum region is then further divided into 
sub-regions. Each sub-region together with its adjacent sub- 
regions as buffer form a subsystem that is simulated by the 
modified neglect of diatomic overlap (MNDO) semi-empirical 
quantum mechanical method (here the MOPAC software had 
been used). The excess charges in the quantum region but 
not in the subsystem being simulated are treated as point 
charges. Their contribution to the subsystem being simulated 
is through Coulomb interaction. 

The coupling of the quantum region and the semi-classical 
region is through the quasi-thermodynamic equilibrium con- 
dition which assumes that the electrochemical potential 
(Fermi level) is a constant along the nanotube, and equal to 
that of the substrate. The densities of the excess charge ("ex- 
cess density" for simplicity) calculated separately in quantum 



3 




(a) (b) (c) 



FIG. 2: Three ending structures of the (5, 5) SWCNT. The gray balls 
stand for carbon atoms, (a) The a- type structure with -BH as ad- 
sorbates. The green balls stand for boron atoms and white balls for 
hydrogen atoms, (b) The b-type structure with -NH as adsorbates. 
The blue balls and the white balls denote nitrogen atoms and hydro- 
gen atoms, respectively, (c) The c-type structure with oxygen ions 
as adsorbates. The red balls stand for oxygen atoms. The Z-axis in 
each figure is parallel to the axis of the tube. The X-axis is the ra- 
dial axis on the plane perpendicular to the tube axis and through a 
hydrogen atom (for FIG. |2(a)| and |2(b)| ) or through the middle point 
of two adjacent oxygen atoms (for FIG. |2(c)l ). The angle a is on the 
X-Z plane. 



and semi-classical regions should coincide at some overlap 
place of two regions. Here we required that the excess densi- 
ties coincide at the position of 900 nm from the substrate. The 
self-consistent excess density of the entire CNT is achieved 
through iterations that contain a small loop and a big loop. In 
the small loop, the sub-regions of the quantum region are sim- 
ulated one by one, and repeated until that a converged elec- 
tron density of the quantum region is obtained. In the big loop 
the quantum region and the semi-classical region are simu- 
lated alternatively until the self-consistent charge distribution 
is achieved. This process is greatly accelerated by the obser- 
vation that the charge density in the region far from the apex, 
to a good approximation, is a linear function of the distance 
from the substrate. This linear behavior is a consequence of 
the constant DOS in the classical region. Bill 

When the SWCNT is mounted on the metallic substrate, 
Schottky junction will be formed at the back contact in prin- 
ciple. However, as the SWCNT here is very long, the effect of 
this junction could be ignored. Equivalently, we have assumed 
that the Fermi level of the isolated SWCNT aligns with that 
of the metal. The boundary condition of the metal surface is 
guaranteed by the image charges of the excess charges of the 
SWCNT. The total electrostatic potential is the superposition 
of the Coulomb potential created by excess charges and their 
images in the substrate, as well as the applied macroscopic 
field. 

The apex structures with -BH, -NH, and -O terminations 
are shown in FIG.[2l they are referred to as a- type, b-type, and 
c-type structures respectively. The five green balls in FIG. 
|2(a)| for the a- type structure stand for the boron atoms. Each 
boron atom shares two bonds with two carbon atoms. The 
third bond of boron atom is saturated by a hydrogen atom (the 
white ball). Similarly, the b-type structure is shown in FIG. 
|2(b)[ with the blue balls for nitrogen atoms and white balls 
for hydrogen atoms. The five oxygen atoms are represented 
by red balls in FIG. |2(c)| for the c-type structure. All ending 
structures are relaxed with the MOPAC software. The coor- 



TABLE II: The coordinates of the adsorbates and the carbons of the 
top layer. Where Z stands for the axial coordinate with the origin on 
the substrate surface, X and Y denote two orthogonal coordinates of 
the plane perpendicular to the tube axis, with the origin at the axis. 



Structures/atoms 


Coordinate (angstrom) 








Z 


X 


Y 


-BH 


H 
B 

C 


10088.04 
10086.92 
10086.04 


2.35 
2.10 
0.94 


2.85 
2.53 
3.31 


-NH 


H 

N 
C 


10087.60 
10086.87 
10086.04 


2.35 
1.91 
0.94 


2.83 
2.31 
3.31 


-0 




c 


10086.83 
10086.04 


2.05 
0.94 


2.46 
3.31 


TABLE III: The binding energy of adsorbates. 




ions 


-BH -NH 


-0 


-H 


binding energy (eV) 


17.8 18.3 


3.99 


8.63 



dinates of adsorbent atoms and the first layer of carbons are 
given in TABLE [III 

The binding energy Eu n d is defined as 

Ebind ~ E se p — E] J0nc i, (1) 

where E sep is the total energy of the system in which the ad- 
sorbates are separated far away from each other and from the 
CNT, and E\j 0n( i is the total energy of the system in which the 
adsorbates are bonding with the carbons of the top layer. In 
the calculation of binding energies, we only consider the last 
ten layers of carbons at the tip. The stability of the ending 
structure is described by the binding energy per adsorbent di- 
atom of the a-type/b-type structure or per adsorbent atom for 
the c-type structure. The results are given in TABLE Hill The 
a-type and b-type ending structures are more stable than the c- 
type structure and the structure in which every dangling bond 
of carbon is saturated by a hydrogen atom. 

III. ELECTRONIC STRUCTURES IN THE ABSENCE OF 
APPLIED FIELD 

The FE characteristic of each individual SWCNT should 
be relied on its intrinsic electronic structure. Therefore in this 
section we focus on the electronic structures of the SWCNT 
in the absence of external macroscopic field. 

FIG. 3(a)/(b)/(c) shows the electron density in the in- 
tersection plane, on which the axis of SWCNT and one 
boron/nitrogen/oxygen atom are located. The arrows mark the 
positions of the adsorbent atoms. In the a-type structure (FIG. 
|3(a)| ), electrons obviously transfer from the boron atoms to the 
hydrogen atoms. From FIG. |3(b)| one can see that electrons 
in the b-type structure are concentrated at the nitrogen atoms. 
The electrons are strongly concentrated at the oxygen atoms 
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TABLE IV: The distribution of the excess charges. 



atom group 


excess charges (e) 


CBH 
CNH 
CO 


-0.082(C), 0.276(B), -0.039(H) 
0.086(C),-0.248(N), 0.150(H) 
0.182(C),-0.302(O) 


TABLE V: The dipoles in the first and the second top layers. 


Structures 


dipole (e-nm) 
atom of 1st layer atom of 2nd layer 


-BH 
-NH 
-0 


-0.0095 0.0135 
0.0158 -0.0122 
-0.0154 -0.0051 



in the c-type structure (FIG. |3(c)) . The electron distribution in 
the c-type structure can be seen more clearly in FIG. ffl where 
the distribution of excess electrons is shown along the wall 
of the SWCNT. The excess charge associated to each atom 
of the top three(two) layers in the a-type/b-type(c-type) struc- 
ture has been calculated by the Mulliken population analysis 
method. The results are given in TABLE HVl Since the boron 
has smaller electronegativity in comparing with the carbon 
and the hydrogen, electrons tend to transfer from the boron to 
both the carbon and the hydrogen in the a-type structure and 
lead to a dipole with its negative charge outward. In the c-type 
structure, we find a dipole of the same direction as the dipole 
in the a-type structure with a bigger value. This is expected 
as Xs of the oxygen is much larger than that of the carbon. A 
reversed dipole, i.e., with its positive charge outward, arises 
in the b-type structure due to the large electronegativity of the 
nitrogen. The dipole values of the top three (two) layers of the 
a-type and b-type (c-type) structures are given in TABLE IVl 



IV. FIELD-DEPENDING BARRIER 

Now let us investigate the response of the SWCNT to the 
macroscopic fields (F app i). The general picture is as fol- 
lows. The applied field F app i drives electrons to the tip of the 
SWCNT until the equilibrium is set up. We assume that the 
emission current is weak and can be ignored in the calcula- 
tion of the electron density of the SWCNT. The redistribution 
of electrons in F app \ has two consequences. First, the field 
is more or less shielded in the body of the tube. This leads 
to the field enhancement at the apex. However, as the wall 
of the SWCNT has only one layer of atoms, the shielding 
would be not complete, especially in the apex region where 
the field is strong. There would be a field penetration in the 
apex region, so the field enhancement factor is smaller than 
that predicted by the classical theory for a metal rod. Second, 
there are excess charges accumulating along the tube, espe- 
cially in the apex region. In the quasi-equilibrium assump- 
tion, the charge accumulation is possible only if the neutrality 
energy level of the SWCNT bends down. This happens when 
the applied field lowers the energies of local orbitals so that 




(c) 



FIG. 3: The distribution of electrons on the intersection plane cross- 
ing one adsorbate atom and the axis of the SWCNT. FIG. |3(a)||3(b)[ 
and |3(c)| are corresponding to the a-type, b-type, and c-type struc- 
tures respectively. 



the orbitals with their energies lying between the Fermi level 
of the substrate and the neutrality energy level can accommo- 
date electrons and contribute to excess charges. 

The superposition of Coulomb potential of all charges in the 
tube and their images in the substrate, together with the poten- 
tial of F app u determines the AVB for the electron emission. In 
FIG. Oa), we plot the electrostatic potential U(z) of the a-type 
structure for various F app \. The Z axis has its origin at the last 
atom of the SWCNT and is parallel to the direction of tube 
axis. The dependence of AVB on F app i is remarkable. The 
barrier potential U(z) of the three ending structures are com- 
pared in FIG. [3b) for 7.0 V/yum. Comparing to the AVB of 
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FIG. 4: The electron distribution around the atoms of last two layers 
of the SWCNT terminated with oxygen atoms. 
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Z (angstrom) 
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10090 10095 
Z (angstrom) 

(b) 



FIG. 5: The AVB of the SWCNT. (a) U(z) of the a-type structure 
with F app i from 6.0 Y/fim to 15.0 V/ji/m, by step of 1.0 V/fim. The 
origin is at one of the last nucleus of the SWCNT. The lower barrier 
is corresponding to larger F app i. (b) U(z) of the three ending struc- 
tures in the applied field of 7.0 V/yt/m. The solid, dashed, and dotted 
curves are corresponding to the a-type, b-type, and c-type structures, 
respectively. Here the origin is on the surface of the substrate. 



the a-type structure, the AVB of the b-type structure is lower 
but thicker. This partly results from the difference of the apex 
dipoles (see TABLE V). Suppose that the applied field would 
not change significantly the intrinsic electronic structure of the 
SWCNT. Then in the a-type (b-type) structure, the atoms of 
the first layer have negative (positive) dipoles, and the atoms 
of the second layer have positive (negative) dipoles. The po- 
tential close to the apex is mostly affected by the atoms of the 
first layer and is enhanced (suppressed). Apart from the apex 
in some distance (7 angstroms), the dipoles of the first layer 
would be screened by those of the second layer, the potential 
of the a-type structure could not be distinguished from that of 
the b-type. For the c-type structure, the dipoles of the first 
two layers are negative, the net negative dipole and the large 
negative excess charge at the apex make the AVB high and 
thick. 

The field enhancement factor as a feature of the AVB is 
highly concerned in the FE applications. For the SWCNT, 
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FIG. 6: The field enhancement factor versus the inverse of applied 
field. The squares, circles, and triangles are data of the simulations 
of the a-type, b-type, and c-type structures, respectively. 



however, it is difficult to define a quantity that corresponds 
exactly to the field enhancement factor as usually defined for 
planar metallic emitters. Denote the maximum local field by 
F apex . Here the field enhancement factor y q is defined as 
F apex/ F app i- The variation of this field enhancement factor 
with the value of l/F app i is shown in FIG. Note that y q 
of these structures do not change significantly with applied 
fields, but are obviously different from each other. It reflects 
the important influence of the adsorbates on FE. 



V. EMISSION CURRENTS AND EMISSION PATHS 

According to the quantum mechanics, electrons may emit 
to the vacuum along any path with non-zero probability. In 
the semi-classical approximation, only the path of the least 
(maximum) action is considered. This path will be referred 
to as the most probable path (MPP). Since the a-type and b- 
type ending structures have the thinnest AVB in front of the 
hydrogen atoms, MPP of these structures should start from 
the hydrogen atom. The thinnest AVB of the c-type structure 
starts from the top carbon atoms and goes outward along a 
path through the middle point of two adjacent oxygen atoms. 
For simplicity, we assumed that the electrons go through the 
barrier in straight lines. The path angles (a) of three ending 
structures are presented schematically in FIG. [2 Along the 
MPP, the transmission coefficient (D) can be estimated by the 
WKB approximation 



Z) = exp[-^J V2mV(z)dz], (2) 

where V(z) is the electron energy potential related to the Fermi 
level, and the integral is over the classical forbidden region, 
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FIG. 7: D versus path angle for the a-type structure in 12.0 Y/fim. 
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FIG. 8: The angle of the most probable path versus the applied field. 
The squares, circles, and triangles are the data corresponding to the 
a-type, b-type, and c-type structures respectively. 



i.e., where V(z) > 0. We has assumed that the electrons pos- 
sess the Fermi level. If the path is not the MPP, the value 
of D by Eq. [2] should be smaller. For instance, In FIG. [7J 
D is plotted against the path angle for the a-type structure in 
12.0 V/yum. The angle of the MPP is determined by the maxi- 
mum of D. Since the AVB is field-dependent, the MPP would 
change its orientation with the applied fields. The angles of 
the MPP versus the applied fields are plotted in FIG. [8j The 
squares, circles, and triangles correspond to the a-type, b-type, 
and c-type structures respectively. One sees that the angle of 
the MPP decreases as the applied field increases. This phe- 
nomenon would be related to the fact that higher applied field 
induces more excess electrons at the apex. 

With D in hand, the emission current is estimated by 



FIG. 9: The currents versus applied field. The squares, the circles, 
and the triangles are the currents corresponding to the a-type, b-type, 
and c-type structures, respectively. The inset is the Fowler-Nordheim 
plot where l/F app i is the horizontal axis and log(I/F^ p/ ) the vertical 
axis. 



(3) 



where q exc are the excess electrons around the atoms from 
which the MPPs start, and v is the collision frequency that can 
be estimated from the average kinetic energy of 7r* electrons 
as Ek(n*)/h. Another way to estimate the collision frequency 
is to use the uncertainty relation 



h 



32n 2 m < Ar 2 > ' 



(4) 



where < Ar 2 > is the uncertainty of the radial coordinate. Its 
numerical value can be estimated from the density of excess 
electrons (FIG. [3]). The collision frequency estimated by two 
methods has the same order of 10 14 Hz. 

The diagram of the emission current versus applied field is 
presented in FIG. [3 where squares, circles, and triangles are 
for the a-type, b-type, and c-type structures, respectively. For 
each structure, the emission current is the summation of the 
currents along five MPPs, each MPP is as defined in FIG. [2 

To confirm that the emission is dominated by the MPPs, 
we have also calculated the emission current along the path 
starting from a charge center in between two nearby MPPs. 
The results for F app i = 12.0 Y/jjm are compared in TABLE 
VI. The contribution of the MPP is one order or so larger than 
that of the path not a MPP. 



VI. CONCLUSIONS AND DISCUSSIONS 

A one micrometer long (5, 5) single- walled carbon nan- 
otube (SWCNT) has been simulated by the multi-scale 
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TABLE VI: The emission currents (A) from the MPP and from the 
path in between two MPPs. 



Structures 


-BH 


-NH 


-0 


MPP 


1.4E-7 


1.5E-7 


1.3E-8 


current 








path not a MPP 


8.7E-9 


1.4E-8 


2.1E-9 



method, with the interest focused at the effects of ending struc- 
tures. The apex- vacuum barriers of the SWCNTs ended by 
-BH, -NH, and -O respectively have been obtained. The local 
field enhancement factor is different from one ending struc- 
ture to another; all are much smaller than the prediction of 
the classical model. The lowering of the apex- vacuum barrier 
by the applied field, besides the local field enhancement, as 
the essential mechanism is confirmed to be responsible for the 
low turn-on field of SWCNTs. 

According to our simulations, the SWCNTs terminated 
with diatoms whose outer atoms have smaller electronegativ- 
ity would be superior in field emission and have more stable 
in structure. The turn-on field ( 77.5 V/yum) of the SWCNT 
terminated with -BH or -NH is smaller than that terminated 
by oxygen. In the same applied field, the current of oxygen 
terminated SWCNT is an order or more weaker than those 
terminated by -BH and -NH. The simulation suggests that the 
binding energy of the diatom decorated SWCNTs is about 4.5 
times the binding energy of the oxygen saturated SWCNT and 
more than double of the hydrogen saturated one. The elec- 
tronegativity of the adsorbates may partly explain the feature 
of the apex- vacuum barrier and give a general hint for the 
structure optimization. The orientation of the most probable 
emission path depends on the applied field significantly. This 
observation would provide a mean for extracting the atomic 
features of the apex. 

One would have noticed that the best turn-on field we esti- 
mated for the -BH ending structure is still too large than most 
experimental results that is about 7-to-5 V/yum. [13211 There 
would be a number of reasons, besides the uncontrollable 
error of the simulation. First, the emission current could 
be larger if one takes into account the emission from the 
paths other than the most probable paths we have considered. 
Secondly, in our calculation, we have ignored the exchange 
and correlation effect (i.e., the image force) on the emitting 
electron. The prediction would be improved by including 
this effect. However, to our knowledge, it is still difficult 
to include the exchange and correlation effect in a quantum 
mechanical simulation of emission current of SWCNT. ll33ll 
Last but not least, most, if not all, CNTs used in the FE 
experiments are multi-walled CNTs and the chirality is not 
controllable. The experimental observed emission current 
would be higher since the multi-walled CNTs would have 
larger field enhancement factor as it would resemble more a 
classical metal rod. The above argument implies the existence 
of discrepancy between our simulational result and the 
experimental observations. With optimism, our simulation 
would have revealed some virtual elements involved in the 
field emission of SWCNT with the atomic decoration, at least 



qualitatively. 
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